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for uniform convergence on compact subsets of C \ [— 1 , 1] . The present pa- 
per is contains a generalization of this fact in the case of families of functions 
satisfying finite recurrence relations with varying functional coefficients. As 
an application we describe the asymptotic ratio for two classes of biorthogonal 
polynomials introduced by Ismail and Masson. The main technical tool we use 
is a parametric generalization of the classical Poincare- Perron theorem on sta- 
bilizing recurrence relations where we assume that the varying coefficients of 
a recurrence depend on auxiliary parameters and converge uniformly in these 
parameters to their limiting values. 



The aim of this paper is to discuss some aspects of the following question. 

Problem. Given a positive integer k > 2 let {<^>j jn (x)}, 1 < i < k, n € Z, be a se- 
quence of fc-tuples of complex-valued functions of a (multi-)variable x = (x\, Xd) 
defined in some domain C R d . Describe the asymptotics when n — * oo of the 
ratio % t 7ff for a family of complex- valued functions {/n(x) | n € Z + } satisfying 



A relation of the form (jl.ip is called a recurrence of length k + 1 or a (k+ 1) -term 
recurrence relation. 

When both the coefficients 0^ n (z) and the functions / n (z) are complex analytic 
one can associate to each meromorphic ratio ^ + ^" > the following useful complex- 
valued distribution v n , see p. 249]. (Since we will only use this distribution in 
one-dimensional case we restrict ourselves to C. For the multi-dimensional case 
consult e.g. [25]-) 

Definition 1. Given a meromorphic function g in some open set f2 C C we con- 
struct its (complex-valued) residue distribution v g as follows. Let {z m \ m £ N} be 
the (finite or infinite) set of all the poles of g in f2. Assume that the Laurent ex- 
pansion of g at z m has the form g(z) = Y2-oo<l<l (z^- ' y ■ Then the distribution 
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u g is given by 

». = E ( E tt^^^^I ' (L2) 

m>l yi</<i m v y 

where S Zm is the Dirac mass at z m . The above sum is meaningful as a distribution 
on O since it is locally finite there. 

Remark 1 . The distribution v g is a complex- valued measure if and only if g has all 
simple poles, see [3] p. 250]. If the latter holds then in the notation of Definition Q] 
the value of this complex measure at z m equals T m ±, i.e., the residue of g at z m . 

Definition 2. Given an integrable complex- valued distribution p in C we define 
its Cauchy transform C p (z) as 




The distribution p itself is then restored from its Cauchy transform C p (z) by 

1 dC p 

P = — 

7T az 

where C p (z) and its ^-derivative are understood as distributions. 

With these definitions we immediately see that any meromorphic function g 
defined in the whole C is the Cauchy transform of its residue distribution v g if the 
condition f c dis g (t;) < oo holds, see p. 261]. 

Definition 3. Given a family {/ n (x)}, x = (xi, ■»,x ( i) of complex- valued functions 
defined in some open set O C l d one calls ^(x) = lim ri _, 00 ^rff the asymptotic 
ratio of the family, provided that ^(x) exists in some open subset of fl. If {/„(z)} 
consists of functions which are analytic in Q and v n denotes the residue distribution 
of the meromorphic function in Q then the limit v = linij^oo v n (if it exists 

in the sense of weak convergence) is called the asymptotic ratio distribution of the 
family {/„(z)}. 

We can now formulate the above problem in a more precise manner: 

Problem specified. Given a family {/„(x)} of functions satisfying a finite re- 
currence relation of the form (jl.ip calculate the asymptotic ratio 'I'(x) and, if 
applicable, the asymptotic ratio distribution v of the family. 

The answers to both parts of this problem are well-known in the case of orthog- 
onal polynomials. (Notice that any family of orthogonal polynomials satisfies a 
3-term recurrence relation of the type (jl.ip with, in general, variable coefficients.) 
Indeed, Szego's fundamental results on orthogonal polynomials assert the following, 
see [21] . 

Theorem 1. Let w(x) be a nonnegative function that satisfies ''^a^^ > ~ °°- 
Then for the sequence of monic polynomials {p n (z) \ degp n (z) = n} which are 
orthogonal on [—1,1] with respect to the weight w(x) one has (in the notation of 
Definitional 

(i) = z+V^l f rzeC\ [-1, 1], 

(ii) v{x) = 2VT ^ dx for x S [-1, 1]. 

The probability measure v{x) in (ii) is supported on [—1,1] and the branch of 
y/ z 2 — 1 in (i) is chosen such that lim 2 ^ oc ^ Z ' ; T 1 = 1. 
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Remark 2. Theorem [T] is actually true for a much broader class of orthogonal 
polynomials, namely, those whose weight w(x) satisfies the requirement w'(x) > 
almost everywhere on the interval [—1,1], see [21]. There are also further general- 
izations of this statement due to P. Nevai, G. L. Lopez and other authors. 

Let us recall some basic facts about linear recurrence relations. Consider first a 
usual recurrence relation of length k + 1 with constant coefficients 

Un+l = Oi\U n + a 2 U n -l + . . . + QtfcUn-fc+1, (1-3) 

with ak 7^ 0. 

Definition 4. The left-hand side of the equation 

t k - {a x t k ~ l + a 2 t k - 2 + . . . + Qk ) = (1.4) 

is called the characteristic polynomial of recurrence l|1.3|) . Denote the roots of (|1.4|) 
by T\, ... ,Tk and call them the spectral numbers of the recurrence. 

Definition 5. Recurrence relation (|1.3|) and its characteristic polynomial Ij 1.4)1 are 
called maxmod- generic if there exists a unique and simple spectral number r maa; 
of this recurrence satisfying |r mQ;r | = maxi<i<fe \n\. Otherwise (|1.3|) and (|1.4p are 
called maxmod-nongeneric. The number r maa; will be referred to as the leading 
spectral number of (|1.3p or the leading root of (|1.4p . 

Let Po/ fc = {t k + ait^ 1 + . . . + a k \ a t E C, 1 < i < k] denote the set of all 
monic polynomials of degree k with complex coefficients, which we identify with 
the fc-dimensional space of complex fc-vectors (ai, . . . , a&). 

Definition 6. The subset C Polk consisting of all maxmod-nongeneric polyno- 
mials is called the standard maxmod- discriminant. For any family 

r(f,ai,. ..,a q ) = {t k + oi(ai, . . . ,a q )t k ~ l + . . . + a fe (ai, . . . ,a q )} 

of monic degree k polynomials in t we define the induced maxmod- discriminant Sr 
to be the set of all parameter values (ai, . . . , a q ) E C q for which the corresponding 
polynomial in t is maxmod-nongeneric. 

Similar discriminants were considered in [T] and )14j . Some local properties of 
Sfe can also be derived from [5] and [17] . 

Example. For k = 2 the maxmod-discriminant S2 C P0I2 is the real hypersurface 
consisting of all solutions to ea\ — 4a 2 = 0, where e E [1, 00) is a real-valued 
parameter. More information on will be given in §6. 

The following simple theorem can be found in e.g. [52 Ch. 4]. 

Theorem 2. Let k E N and consider a k-tuple (ax, . . . ,ctk) of complex numbers 
with a k 7^ 0. For any function u : Z + — > C the following conditions are equivalent: 

(i) J2 n >Q u n tn = gfi}> where Q(t) = 1 - (otit + a 2 t 2 + . . . + a k t k ) and P(t) is 
a polynomial in t whose degree is smaller than k. 

(ii) For all n > k — 1 the function u n satisfies the (k + \)-term recurrence 
relation (|1.3|) . 

(iii) For all n > one has 

r 

u n = Y,Pi(n)r?, (1.5) 
»=i 

where t\ , . . . , r r are i/ie distinct spectral numbers of (|1.3j) wiit/i multiplicities 
mi, . . . ,m r , respectively, and Pi(n) is a polynomial in the variable n of 
degree at most rrii — 1 /or 1 < i < r. 
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Definition 7. An initial fc-tuple of complex numbers (uq,ui, . . . ,Uk-i) is called 
fast growing with respect to a given maxmod-generic recurrence of the form (|1.3p if 
the coefficient n max of r™ lax in (|1.5|) is nonvanishing, that is, u n — ^max^max + • ■ ■ 
with K max 7^ 0. (Note that by Definition [5] the leading spectral number r max 
of any maxmod-generic recurrence has multiplicity one.) Otherwise the fc-tuple 
(uq, mi, ... , u/s-i) is said to be slow growing. One easily check that the set of all 
slowly growing initial A:-tuples is a (complex) hyperplane in its complement 
being the set of all fast growing fc-tuples, see Lemma [31 The latter hyperplane of 
slow growing fc-tuples can be easily found explicitly using linear algebra. 

A famous and frequently used generalization of Theorem [5] in the case of vari- 
able coefficients was obtained by H. Poincare in 1885, [20] and later extended by 
O. Perron, |19j . It originally uses a continuous variable x, and its discrete version 
follows by restricting x to the set of integers. 

Theorem 3 (Poincare-Perron). If the coefficients Qi(x), i = 1, of a linear 
homogeneous difference equation 

f(x + k) + Qi{x)f(x + k-l) + Q 2 (x)f(x + k - 2) + ... + Qk(x)f[x) = (1.6) 

have limits linx^oo Qi{x) = ctj, i = 1, k and if the roots of the characteristic 
equation \ k + ai\ k ~ 1 + ... + = have distinct absolute values then the limit of 
the ratio ^fe^p f 0T x — > oo of any solution f(x) of (|1.6p equals one of the roots of 
the characteristic equation. 

One can also show that the set of solutions for which lim^^oo ^fe^ — ^max 
where X max is the (unique) root with the maximal absolute value is a (complex) 
hyperplane in the space of all solutions, see Lemma [H But contrary to the case of 
constant coefficients, there seems to be no easy way to determine this hyperplane 
explicitly. 

A number of generalization of the Poincare-Perron theorem can be found in 
the literature, see e.g. [T^], [Hj, [T5], most recent of those being apparently the 
next result obtained by the 2nd author in connection with the present project, [H 
Theorem 1.2]. Let Mk(F) and GLk(F) denote the spaces of all and resp. of all 
invertible (fc x fc)-matrices over a field F. 

Theorem 4. Let {T n }neN be a sequence of regular matrices in GLfc(C) converging 
to some (possibly singular) matrix Ot/Te Mfe(C). Assume furthermore that T 
has a positive spectral radius p(T) and that the circle {z € C | \z\ = yo(T)} contains 
exactly one eigenvalue X of T which is a simple root of its characteristic equation. 
Let u denote an eigenvector of T corresponding to \, i.e., Tu — Au, O^ue C fc . 
Then the complex line spanned by the product T n T n ^\- ■ - T\ € Mfc(C) converges 
to the complex line spanned by uw* € Mfc(C) for some fixed vector 0^w£C l . 
Hence for any vector Xo G C fe such that w*xo 7^ the complex line in C k spanned 
by T„T„_i • • -TiXo converges as n — > 00 to the complex line in C k spanned by u. 

Our first result is the extension of the latter result (and therefore of the original 
Poincare-Perron theorem) to a sequence of invertible matrices depending continu- 
ously on auxiliary parameters. 

Theorem 5. Let {T„(x)}„ e pj be a sequence of families of matrices in GLfc(C) 
depending continuously on x G T> C R d . Assume that this sequence converges 
uniformly, on any compact set in T>, to a matrix T(x). Suppose furthermore 
that for each x G T> , T(x) has exactly one simple eigenvalue A(x) of the maxi- 
mal modulus with the corresponding eigenvector u(x) G P fc_1 . Then the product 
T„(x)T„_i(x) . . . T2(x)Ti(x), viewed as an automorphism of¥ k ~ 1 , converges to the 
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transformation u(x)w t (x) on P fe 1 \H(x), where w(x) G P fe 1 is continuous in T> 
and -ff(x) is the hyperplane given by w'(x)v = 0. 

We can now formulate our further results which are basically corollaries of The- 
orems S] and [5] 

Notation 1. Let k E N and assume that 

{$„ := ($i,„(x),$ 2 ,n(x),. .. ,$fc,„(x)) | n E Z+} , x= (x u ...,x d ) 

is a family of fc-tuples of complex- valued functions defined on some domain SI C K d 
such that <!>„ converges to a fixed fc-tuple of functions ^$i(x), $2(x), . . . , $fc(x)^ 

pointwise in SI. Denote by IT = (/o(x), . . . , /fc_i(x)) some initial fc-tuple of func- 
tions defined on SI and let {/„(x) | n E Z + } be the family of functions that satisfy 
the recurrence relation 

fe 

/ n (x) = $i,n-l(x)/„-i(x), z G O, (1.7) 
i=l 

for all n > A; and coincide with IT = (/o, . . . , /fc-i) for < n < — 1. 

Theorems 2] and [S] imply the following. 

Theorem 6. In the above notation there exists a unique subset C f2\S$ which 
is minimal with respect to inclusion and such that the following holds: 

(i) For any x G SI \ (S$ U £$) one has 

,. /n+l(x) . . 

where VP; (x) is £/ie leading root of the asymptotic symbol equation 

* fe ^$ 1 * fc - 1 + $ 2 * fe - 2 + ... + $ fe (1.8) 



and E$ denotes the induced maxmod- discriminant of (|1.8j) considered as a 
family of monic polynomials in the variable 'J (cf. Definition^). 
(ii) If $ n consists of functions continuous in SI and <&„ =4 ($i, $2, ■ ■ ■ , &k) in 
SI then 

^* ( ( x ) *n fi\(H $ US*) 1 

JnyX-) 

where ^ stands as before for uniform convergence on compact subsets. 

Remark 3. Notice that since we assumed the uniform convergence in (ii) above 
then each $i is continuous in Q. Moreover, since a uniformly converging sequence 
of analytic functions necessarily converges to an analytic function one has that (ii) 
holds in the analytic category as well. 

Definition 8. The set E$ introduced in Theorem [6] is called the set of slowly 
growing initial conditions (cf. Definition [7]) . If the functional coefficients in (|1.7[) 
are fixed (i.e., independent of n) then E$ is exactly the set of all points p G SI such 
that the initial fc-tuple (fo(p), ■ ■ ■ , fk-i(p)) is slowly growing with respect to the 
recurrence (|1.7p evaluated at p, see Definition 6. 

Remark 4. The exact description of in the case of varying coefficients can be 
found in the proof of Theorem [6l see §2. Note that the convergence of —^l- to the 
leading root of (|1.8p might actually occur at some points lying in 5$ where the 
leading root is still unique but multiple, see Lemma[3] (iii)-(iv). 
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Notation 2. Using Notation 1 let V — fl x INk denote the Cartesian product of 
n and the complex fc-dimensional linear space INk of all initial fc-tuples of complex 
numbers. Any initial fc-tuple IT — (/o(x), . . . , /fc_i(x)J of respectively, analytic 
functions in can be considered as an analytic section of V viewed as a trivial 
vector bundle over fl. Denote by V the restriction of the bundle V to the subset 
fl \ S$ of the base. 

Definition 9. Let A$ denote the subset of V consisting of all pairs (p, ink) such 
that the asymptotic ratio of the recurrence (jl.7p evaluated at p with initial fc-tuple 
ink G INk does not coincide with the leading root of ()1.8p at p. The set A$ is 
called the subvariety of slow growth. 

Proposition 1. In the above notation the subvariety A$ of slow growth is an 
analytic subbundle of complex hyperplanes in V provided that the family of ' k -tuples 

{$„ = ($ lin (x),$ 2 ,n(x),...,$fe ) n(x)) | U G 1+ } 

and its limit k-tuple ^$i(x), <£>2(x), . . . , $fc(x)^ consist of analytic functions that 
satisfy <f n =t ($i,$ 2 , . . . , $ fc ) in O. 

Definition 10. Let £ : _E — > B be a (complex) vector bundle over a base i? 
with linear fiber L. A smooth (respectively, analytic) section s : B — > _E is called 
transversal to a smooth (respectively, analytic) submanifold 7i C -E if at each point 
p of the intersection sDH the sum of the tangent spaces at p to s and 7i , respectively, 
coincides with the tangent space at p to the ambient space E. 

Remark 5. Recall that a subset of a given set is called massive if it is the intersection 
of at most countably many open dense subsets. Thorn's transversality theorem, see 
e.g. [5], implies that the set of all smooth sections of £ : E — * B which are transversal 
to a given smooth subbundle H C E is a massive subset of the set of all sections. 
The same holds in the analytic category provided that the fibration is trivial. 

Corollary 1. In the notation of Theorem^ if an initial k-tuple IT of analytic 
functions in O is transversal to the subvariety of slow growth A$ then the set E$ 
of slowly growing initial conditions is either empty or an analytic subvariety in Q 
of real codimension two. 

In particular, in the one-dimensional case one has the following. 

Corollary 2. In the notation of Theorem^ the following holds: If the coefficients 
(3>i >n (z), . . . , $ fe> „(z)) ; n G Z_|_, and the initial k-tuple IT — (fo(z), ■ ■ ■ , fk-i(z)) 
of recurrence (|1.7[) are complex analytic functions in an open set f2 C C then either 
S$ = fl or S$ is a union of real analytic curves while is either the whole fl or 
consists of isolated points, see Fig. 1. 

Recall now the asymptotic ratio distribution v introduced in Definitions [T] and O 

Proposition 2. In the complex analytic category the following holds: 

(i) The set £$ of slowly growing initial conditions is irrelevant to the support 
of v. Thus supp^ C S$, where £$ denotes as before the induced maxmod- 
discriminant. 

(ii) Suppose that there exists a nonisolated point po G S$ such that equation 
(|1.8[) considered at p$ has the property that among its roots with maximal 
absolute value there are at least two with the same maximal multiplicity 

(compare with Definition \11\ below). If the sequence | " I n G 

diverges then the support of v coincides with S$ . 
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Figure 1 . Zeros of polynomials satisfying the 4-term recurrence 
relation p„+i(z) = p n (z) + (z + l)(z - I)p n -i(z) + (z 3 + 10)p„_ 2 {z) 



Explanations to Figure\j\ The nine fat points in these three pictures are the branch- 
ing points of the symbol equation * 3 = * 2 + (z + \){z - /)* + (z 3 + 10). The two 
upper pictures show the zeros of pmiz) and P4s(z) for the initial triple po(z) — 
0, Pi{z) = z 4 — 5/, P2(z) — z. The lower picture shows the zeros of Pqa{z) for the 
same recurrence relation but with the initial triple pa(z) = z s — z 5 + /, pi(z) — 
z — 51, pi{z) = 51 z 2 + z — 10. (Both initial triples were picked randomly.) One 
observes that on all three pictures the zeros form a pattern close to a real analytic 
curve 5$ and a number of isolated points. On the upper pictures there are four iso- 
lated points which practically coincide on both pictures although the polynomials 
themselves are coprime. On the lower picture there are seven isolated points which 
also form a very stable set, as these points hardly change if one takes different poly- 
nomials Pn(z) with the same initial triple. Unfortunately, limitations of the package 
Mathematica™ and of our computer have not allowed us to check how close these 
isolated zeros are to the actual set of slowly growing initial conditions for the 
cases considered. 

Remark 6. The isolated points on Figure [T] have a strong resemblance with the 
spurious poles that were considered in a substantial number of papers on Pade 
approximation, see e.g. [23j and references therein. The study of the exact relation 
between these two objects can be found in jjj. 

This paper is organized as follows. In §2 we prove our generalization of the 
Poincare-Perron theorem and settle the remaining results in §3. In §4 we consider 
concrete examples of 3-term recurrence relations with polynomial coefficients re- 
lated to the theory of biorthogonal polynomials [IT] . In §5 we discuss a number of 
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related topics and open problems. Finally, in appendix we study the topological 
structure of the standard maxmod-discriminant C Polk- 

2. Proving parametric generalization of Poincare-Perron theorem. 

For F = R,C denote by F fe , M fc (F), GL fe (F) the fc-dimensional vector space, the 
algebra of k x k matrices and the subgroup of k x k invertible matrices over the 
field F. Denote by || • || any vector norm on F fe or on Mfc(F). Let || ■ H2 be the 
£2 norm on F fe induced by the standard inner product (x,y) := y*x on F fc and 
denote by || • H2 the induced operator norm on Mfc(F). For T G M^(C) denote by 
p(T) the spectral radius of T, i.e. the maximal modulus of the eigenvalues of T. 
As in \7\, the spaces PF fe ,PM fe (F),PGL fc (F) arc obtained by identifying the orbits 
of the action of F* :— F\{0} on the nonzero elements of the corresponding sets by 
multiplication. Then PM fc , FM k (R) and PC fc , PM fe (C) are compact real and complex 
manifolds respectively. To keep our notation standard we set P fe_1 = PC fc . For 
x 6 F fc \{0},T G M fe (F)\{0} denote by x,f the induced elements in PF fc ,PM fe (F) 
respectively. Furthermore, T can be viewed as a map from PF fc \kerT to PF fe . 
Denote by d(-, •) : P fc_1 x P fe_1 — > [0, 00) the Fubini-Study metric on P fc_1 , see e.g. 

Consider an iteration scheme 

x„ := T n x„_!, x G F fc , T n G M k (¥), n G N. (2.1) 

This system is called convergent if x ra , n G N is a convergent sequence for each xo G 
F fc . This is equivalent to the convergence of the infinite product ...T n T n ^i...T2Ti, 
which is defined as the limit of T n T n -\...TiT\ as n — ► 00. For the stationary case 
T n = T, n G N the necessary and sufficient conditions for convergency are well 
known. Namely, (i) the spectral radius p{T) can not exceed 1; (ii) if p(T) = 1, 
then 1 is an eigenvalue of T and all its Jordan blocks have size 1; (iii) all other 
eigenvalues A of T different from 1 satisfy |A| < 1. 

In some instances, as e.g. for problems related to the Lyapunov exponents in 
dynamical systems, one is interested if the line spanned by the vector Xj converges 
for all Xo ^ in some homogeneous open Zariski set in F fc . If this condition holds 
we call (12.11) projectively convergent. 

For the stationary case ^ T n = T G A/fe(C) one can easily check that (|2.1|) is 
projectively convergent if and only if among all the eigenvalues A of T satisfying 
|A| = p(T), there is exactly one eigenvalue Ao which has Jordan blocks of the 
maximal size. 

The main point of Theorem \E\ is to show that w(x) is continuous in T>. We 
use some technique developed in [7], in particular, the arguments of the proof of 
TheoremlU [71 Theorem 1.2, page 258 ], to obtain the continuity of w(x). Namely, 
let A(x) is the unique simple eigenvalue of T(x) with the maximal modulus. Then 
A(x) is continuous and nonvanishing in T>. To prove Theorem[5j we replace T(x) and 
T n (x) by j^jT(x) and j^)T n (x) for each n, respectively. Thus we can assume 
that for each x G T> the number A = 1 is the simple eigenvalue of T(x) of the 
maximal modulus. Let u(x),v(x), 7^ for each x G T>, be the right and resp. the 
left eigenvectors of T(x) corresponding to 1, i.e. 

T(x)u(x) = u(x), T t (x)v(x) = v(x), u t (x)v(x) = 1, xeD. (2.2) 

Since u(x),v(x) can be chosen continuous in T> it follows that u(x),v(x) G P" _1 
are continuous in T> as well. In what follows we fix a norm || • || on Af^(C) satisfying 
for any T G M fc (C) the condition ||T|| = ||T rt |] . 

We need two additional Lemmas; the next one generalizing the inequalities given 
on pages 254-255 of [7J. 
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Lemma 1. Let {T n (x)} ri gN be a sequence of(kx k) complex- valued invertible matri- 
ces, with continuous entries for xePc Assume that this sequence converges 
uniformly, on any compact set in T>, to a matrix T(x). Suppose furthermore that 
for each x 6 T>, the number A = 1 is a simple eigenvalue o/T'(x), and all other 
eigenvalues o/T(x) lie in the open unit disk |A| < 1. Then for any fixed Xo and for 
each e > t/iere exists 8 — 6(e) > and N — N(e), m — 771(e) S N such that 

\\T n+m {x.)T n+m _ x {x) . . . T„+i(x) - u(x )v t (x )|| < e, (2.3) 
/or eac/i n > N and |x — xo| < <5. 

Proof. Since 1 a simple eigenvalue of T(xo) and all other eigenvalues of 
T(xq) lie in the open unit disk, it follows that linim^oo T m (xo) = u(xo)v t (xo). 
Hence, there exists m = m(e) so that ||T m (x ) — u(x )v t (x )|| < |. Since the 
product map (Mk(C)) m — * Mfc(C) sending a m-tuple X\, X m to their product is 
continuous, there exists Si > 0, such that for — T(xo)|| < S% for i = 1, . . . , m the 
inequality HX1X2 . . . X m — T(x ) m || < § holds. Since T„(x) converges uniformly 
to the continuous T(x) on any compact subset of T>, it follows that there exists 
5 = (5(e) > and N — N(e) G N such that ||T n (x) - T(x )|| < | for n > N and 
||x-x || < S. □ 
Next we give the following modification of [3 Lemma 5.1]. 

Lemma 2. Let E G Affc(C) be a matrix of rank one with p(E) > 0, i.e. E = 
vu\ u'v =1. Set O r := {x G PC fc : rf(x, v) < r} such that O r n kcri? = 0. Then 
E : O r — * {v}. T/ien f/iere exists e — e{r) such that the following conditions hold. 
Assume that B G Mfe(C) satisfies \\B — E\\ < e. Then 

(1) d(BZ,v) < § for each x G O r . 

(2) d(B±,By) < 5d(x,y) /or eac/i x,y G O r . 

Proof. Clearly _Ex = v if u'x 7^ 0. Hence £ : O r — > {v}. Recall that for any 
fl G M fc (C)\{0} the transformation B : PC fc \keri? -> PC fc is analytic. Hence there 
exists e > so that for ||5 - E\\ < £0, B ^ 0, and kerS n O r = 0. Thus, B is 
analytic on O r . One has that for a small e G (0, eo) any B satisfying \\B — E\\ < £ 
is a small perturbation of E. Since i?O r — {v} we deduce condition [1] for any 
e G (0, £i) for some < £\ < £q. Since d(Ex, Ey) = for any x, y G O r we deduce 
condition [5] for some e G (0, £1). □ 

Proof of Theorem \S[ Recall that [7, Theorem 1.2] implies that the map 

T n (x)T n _i(x) . . . T2(x)Ti(x), viewed as an automorphism of P fc_1 , converges point- 
wise to the transformation u(x)w'(x) on P k ^ 1 \H (x) . Now fix some xo G T>. Then 
there exists 8\ > such that the eigenvectors u(x),v(x) G C k satisfying ()2.2|) can 
be chosen continuously in |x — xo| < Si. Since v*(x)u(a;) = 1 it follows that we can 
choose 82 G (0, Si) such that v(xo)u(x) ^ for |x— xo | < 61. To show the continuity 
of w(x) in |x — xo| < S for some S3 G (0, ^2) it is enough to show that the continuity 

of the limit of {^^)TJ^)tZ^J!x.) ■ ■ ■ T^{^T\{xj} in |x - x | < S 3 . Since each 
Tj(x) G GL(fe,C) is continuous, it is enough to show the continuity of the limit of 

{v t (x )r n+ Ar(x)T I1+A r_ 1 (x) . . . Tj V +2(x)T w+ i(x)} in |x-x | < ^^One has to show 

that a subsequence of {v t (x )T„ + jv(x)T„ + Ar_i(x) . . . Tjv+2(x)Tjv+i(x)} converges 
uniformly in |x — xo| < ^3 in the Fubini-Study metric d(-, ■) on P* -1 . 

This is done by using Lemmas QJH Namely, let E = v(xo)u t (xo). Choose 
r > 0, £ > so that the conditions [TJ2] of Lemma [5] hold. Next choose N, S so that 
the condition (|2.3p holds. Let £3 = min(<5, 62). Define 

-^j( X ) = ^(j-l)m+JV+l( X ) ' ' • ^im+w( x ) f° r i = 1) 
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Assume that |x — xo| < £3. We claim that the sequence {i?i(x) . . . £?;(x)v(xo)} 
converges uniformly to w(x) G O r in the Fubini-Study metric. First observe that 
(|2.3|) yields that Bj(x)O r C Or, by condition Q] of Lemma [2 Hence y;(x) = 
Bi(x) . . . £?/(x)v(xo) C Or. Condition [T] of Lemma [U yields that 

d{B^j . . . B^j^oj, S^x) . . . s7^x)vR)) < 

< ^d(^o), bQx) . . . fC^W^R)) < j v 

Hence the sequence {yz}i<=N converges uniformly in |x — Xo| < £3 to y(x). Since 
each y;(x) is continuous in |x — x | < 5% it follows that y(x) is continuous in the 
open disk |x — xq| < S3. □ 

3. Proving remaining results 

We start with a simple consequence of Theorem [2 

Notation 3. Let Reck be the fc-dimensional complex linear space consisting of 
all (fc + l)-term recurrence relations with constant coefficients of the form (|1.3[) . 
As before we denote by INk the ^-dimensional complex linear space of all initial 
fc-tuples (uq, ...,Ufc_i). 

Definition 11. A maxmod-nongeneric recurrence relation in Reck with initial fc- 
tuple ink € INk is said to be of dominant type if the following conditions are 
satisfied. Let t±, . . . ,T r , r < fc, denote all distinct spectral numbers with maximal 
absolute value and assume that these have multiplicities mi , . . . , m r , respectively. 
Then there exists a unique index io G {l,...,r} such that nii < TOj for i € 
{1, . . . , r} \ {«o} and the initial fc-tuple ink is fast growing in the sense that the 
degree of the polynomial Pi in (|1.5p corresponding to t; is precisely mj — 1. The 
number Ti is called the dominant spectral number of this recurrence relation. 

Lemma 3. In the above notation the following is true: 

(i) The set of all slowly growing initial k-tuples with respect to a given maxmod- 
generic recurrence relation in Reck is a complex hyperplane SQk in INk- 
The set SQk is called the hyperplane of slow growth. 

(ii) For any maxmod- generic recurrence relation in Reck and any fast growing 
initial k-tuple (uq, . . . , Uk—i) the limit limn^oo exists and coincides 
with the leading spectral number T max , that is, the (unique) root of the 
characteristic equation (|f .4[) with maximal absolute value. 

(hi) Given a maxmod-nongeneric recurrence relation of dominant type in Reck 
the limit linin^oo exists and coincides with the dominant spectral num- 
ber. 

(iv) For any maxmod-nongeneric recurrence relation of nondominant type in 
Reck the set of initial k-tuples for which lirriji^oo exists is a union of 
complex subspaces of INk of positive codimensions. This union is called the 
exceptional variety. 

Proof. In order to prove (i) notice that the coefficient K max in Definition [7] 
is a nontrivial linear combination of the entries of the initial fc-tuple (uq, . . . , Uk-i) 
with coefficients depending on a±, . . . , otk- Therefore, the condition n max = de- 
termines a complex hyperplane SQk in INk- One can easily see that the hyperplane 
of slow growth is the direct sum of all Jordan blocks corresponding to the spectral 
numbers of a given recurrence (|1.3[) other than the leading one. 

The assumptions of part (ii) together with (|1.5p yield u n — KmaxT 1 ^^ + ■ ■ ■ for 
n £ , where the dots stand for the remaining terms in JT3]) corresponding to the 
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spectral numbers whose absolute values are strictly smaller than |r max |. Therefore, 
the quotient has a limit as n — > oo and this limit coincides with T max , as 

required. By definition T max is a root of (|1.4|) . which completes the proof of (ii). 
This last step can alternatively be carried out by dividing both sides of (|1.3p by 
U n -k+i and then letting n — > oo. In view of Definition [TT] the same arguments show 
that the assertion in (iii) is true as well. 

For the proof of part (iv) we proceed as follows. Take any maxmod-nongeneric 
recurrence relation of the form (|1.3|) and let r±, . . . , r r , r < k, be all its distinct 
spectral numbers with maximal absolute value. Thus \Ti\ = \T max \ if and only 
1 < i < r. Choose an initial fc-tuple IT = (v,q, . . . , Uk—i) and denote by Pi, P r 
the polynomials in (|1.5[) corresponding to n, . . . ,r r , respectively, for the sequence 
{u n | n € constructed using the given recurrence with initial fc-tuplc IT as 
above. Assuming, as we may, that our recurrence relation is nontrivial we get 
from (|1 .5|) that ^ if 1 < i < r. We may further assume that the degrees 
d±, . . . , d r of the polynomials Pi, P r , respectively, satisfy d\ > . . . > d r . A direct 
check analogous to the proof of part (ii) shows that if only the polynomial Pi 
is nonvanishing then linin^oo = t\. If r > 2 and s G {2, . . . , r} is such that 
Pi , . . . , P s are all nonvanishing polynomials among P±, . . . ,P r then using again (|1.5|) 
we get 

(\ n+l / \ n+1 
— J +... + P s (n+l) ( — ] +o(l) 
ny \nj 

Un Pi(n) + P 2 (n) (^\ +... + P s (n)(^\ +o(l) 

Since |^-| = 1 and r,- 7^ n, 2 < i < s, it follows that if di = di then the expression 
in the right-hand side has no limit as n — > 00. Therefore, if such a limit exists then 
di > g?2 , which gives us a complex subspace of /TV^ of (positive) codimension equal 
to d\ — g?2 ■ Thus the exceptional variety is a union of complex subspaces of IN^ of 
(in general) different codimcnsions. 

Part (i) of Lemma [3] is actually a special case of Theorem 2.2]. □ 

Let us now present analogs of this lemma for recurrence relations with varying 
coefficients of the form 

U n + 1 = Oii. n U n + a 2 ,„U n -l + . . . + akjnU n -k+l (3.1) 

where ak, n £ C, n > k — 1. 

Lemma 4. Assume that for all i S {1, . . . , k} the limit linin^oo ai^ n =: on exists 
and denote the limiting recurrence relation of (|3.ip by 

u n +i = a>iu n + a 2 u n -i + . . . + a k u n - k +\- (3.2) 

If (|3.2[) is maxmod- generic then limn—^ exists and equals T max on the com- 

plement of a complex hyperplane SQk in IN^, where T ma x denotes the leading spec- 
tral number of (|3.2p . By analogy with the case of a recurrence relation in Reck 
we call the exceptional complex hyperplane SQk the hyperplane of slow growth 
(cf. Lemma\B$. 



12 



J. BORCEA, S. FRIEDLAND, AND B. SHAPIRO 



Proof. The result is an immediate corollary of Theorem U] applied to the 
family of linear operators {T n \ n £ N} whose action on the space INk is given by 



V o 



Ot2,i 



1 







<*3,? 






ot k , 





A 



i 



n £ N. 



/ 



□ 



From Theorem H] we deduce the following corollary (cf. Proposition 2.1]): 

Corollary 3. Under the assumptions of Theorem^ there exists a sequence of real 
numbers {O n } n efi such that 

T n T n ~i • • • Ti 



lim e l 



•711 



uw 



the convergence taking place in (C) . 

Proofs of Theorem [5] and Proposition [1] are now applications of the above results 
together with Theorem [5] 

Proof. [Theorem [5] (i)] Let us first determine £$ and at the same time prove 
the pointwise convergence of to ^i(xi, . . . , Xd) in the complement ft \ 

(S$ U £$). One can view of the set-up of Theorem [H] as the situation when the 
recurrence relation (|3.ip depends on additional parameters (xi, . . . ,Xd)- Thus at 
a given point p £ Q C R d the limit lim^oo ^rrffi does not exist in one of two 
possible cases. In the first case, p lies in the induced maxmod-discriminant S$ 
(cf. Definition [5]) and the initial fc-tuple (fo(p), ■ ■ ■ , fk~i{p)) is sufficiently generic. 
In the case of fixed coefficients one can be more precise and say that the initial fc- 
tuple (fo(p), ■ ■ ■ i /fc-i(p)) avoids the exceptional variety defined in Lemma [3] (iv). 
(Unfortunately, the analog of this statement for recurrence relations with varying 
coefficients is at the moment unknown, see §5.) In the second case, p lies in \ S$ 
and the corresponding initial fc-tuple (fo(p),---,fk-i{p)) lies in the hyperplane 
SGk (p) of slow growth at the given point p, see Lemma |3J The latter set of points 
p is by definition £$ and it is called the set of slowly growing initial conditions. In 
the case of fixed coefficients in (|1.7[> one can again be more precise: one has a fixed 
recurrence relation in Reck at each point p £ f2\S$. The hyperplane of slow growth 
at p is then the sum of all Jordan blocks corresponding to all nonleading spectral 
numbers of the considered recurrence (cf. the proof of Lemma (i))- Therefore, 
in this case £$ is just the set of all points p £ \ S$ described in Definition [71 
For varying coefficients the hyperplane of slow growth at a given point p £ O \ S$ 



is determined by lim,, 



„ie„ t„t„- 



see Corollary [31 Indeed, Theorem [H 



° ||T n T n _i-2U__ 

implies that in this case the hyperplane SGk (p) of slow growth at p consists of all 
vectors x(p) £ C k such that w(p)*x(p) = 0. 

To settle Theorem [S] (ii) we now assume that <!>„ =4 ($i, ( &2, ■ ■ • , ^fe) and all 
considered functions are continuous in fl. Let po £ f2\ (5$US$). From Corollary[3] 
and the arguments of §4 and §5] it follows that one can find a sufficiently small 
neighborhood O po C Q\ (S$ U S$) of po such that for any e > there exists N e £ N 
satisfying 

T n (p) 



ie„{p)_ 



u(p)w*(p) 



< e for p £ O po , n > N t , 



(3.3) 



\\Tn( P )\\ 

where T n (p) := T n (p)T n _i(p) • • • T\[p). Clearly, u(p)w t (p) is a rank one matrix 
and so it has precisely one simple eigenvalue of maximum modulus for all p £ 
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O po . From (|3.3[) we deduce that there exists N 6 N such that the polynomial 
det(A/fc — T n (p)) has precisely one simple eigenvalue X n (p) of maximum modulus 
for all p € O po and n > N. The implicit function theorem implies that X n (j>) is 

continuous in Pn for all n > N, so that < " > is a sequence of matrix- 

U»(p)Jn>JV 

valued functions which are continuous in Pa . On the other hand, it is not difficult 
to see that this sequence converges pointwise in O po to u(p)w*(p) (see e.g. the proof 

of Theorem H given in [6 ). Together with J33J) this implies that ( tv44t 1 

I \ X n(P)\ J n > N 

is a sequence of (continuous) functions that converges pointwise to the constant 
function 1 on O po . It follows that there exists some open subset O pQ C O po with 

Po € O' such that / tt^tt f is a uniformly bounded sequence of functions 
on 0' pg . Thus <| ~ ~ ) N J- is a bounded sequence of continuous matrix- valued 



Kip) 



i>N 

functions which converges pointwise to the function u(p)w'(p) in 0' pQ . Invoking 
Vitali's theorem we conclude that 

^^u(pW(p) in O' po . (3.4) 

□ 

Next we verify Corollaries [1] and [21 

Proof. [Corollary [1] By Thorn's transversality theorem the set of all smooth 
initial fc-tuples which are transversal to a given subvariety A$ of slow growth is 
a massive subset of the set of all possible smooth initial fc-tuples (cf. Remark [5]) . 
Combined with Proposition [T] this implies in particular that for a given initial k- 
tuple IT transversal to A$ one has that IT (considered as a section of V) and A$ 
are smooth submanifolds of V transversal to each other, where V is the restriction 
of the bundle V = x INk to the subset fi \ S$ of the base (cf. Notation [2]). 
This transversality property implies in turn that the intersection IT D A$ is either 
empty or a smooth submanifold of V of real dimension d — 2, where d is the real 
dimension of f2. Notice that the image £(7TnA$) of the projection £ : V — * S1\S$ 
to the base is exactly £$. Actually, £ induces a diffeomorphism between IT n A$ 
and To see this one simply notes that £ maps the section /T diffeomorphically 
to Q, \ 3$ and that the intersection IT n A$ is a smooth submanifold of the section 
IT. The proof of the corresponding result in the analytic category is analogous. 

Proof. [Corollary [2] By Proposition Q] the set A$ is analytic in the analytic 
category. Thus in this case the asymptotic symbol equation (|1.8j) is either maxmod- 
nongeneric everywhere in O C C or on a one-dimensional real analytic subset, which 
proves the statement in part (i) concerning 5$. (Note that in fact S$ is always 
a real analytic set.) Analogously, the set £$ is analytic in the analytic category. 
Therefore, the analytic continuation of is either a complex analytic curve, in 
which case it coincides with f2, or an analytic zero-dimensional subset of Q, i.e., the 
union of isolated points. 

We finally turn to Proposition [2} 

Proof. [Proposition [2] Under the assumptions of the proposition the leading 
root of the asymptotic symbol equation (|1.8p is a well-defined analytic function in 
a sufficiently small neighborhood of S$. Therefore, the residue distribution (|1.2|) 
associated to this leading root vanishes in a neighborhood of £$. Thus the set £$ of 
slowly growing initial conditions can be deleted from the support of the asymptotic 
ratio distribution v, which proves (i). 
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In order to prove that under the nondegeneracy assumptions of (ii) the support 
of v coincides with 5$ we show that in this case the sequence : %7^ diverges for 
almost all p £ S$. Indeed, due to analyticity the first condition of (ii) implies 
that the recurrence relation is of nondominant type almost everywhere in H$. The 
second condition of (ii) then implies that the sequence ^rff diverges almost 
everywhere in S$, which settles part (ii) of the proposition. □ 

4. Application to biorthogonal polynomials 

Let us first recall some definitions from [11J . In op. cit. the authors introduced 
two types of polynomial families related to what they call Ri- and i?//-type con- 
tinued fractions, respectively. These families were later studied in [27] . 

A polynomial family of type Rj is a system of monic polynomials generated by 

Pn+i(z) = (z - c n+ i)p n (z) - X n+ i(z - a n+ x)p n -i(z) 

with p~i(z) = 0, p {z) = 1 and p„(a n+ i)A„ + i ^ for n e Z+. 

A polynomial family of type Ru is a system of monic polynomials generated by 

Pn+i(z) = {z- c n+ i)p n {z) - X n +i(z - a n+ i)(z - b n+1 )p n -x(z) 

with p-i(z) = 0, p (z) = 1 and p n {a n+ i)p n (b n+ i)X n+1 ^ for n e Z+. 
Assuming that 

lim a n — A, lim b n = B, lim c n — C and lim A„ = A, 

n — >oo n — *oo n — >oo n — >oo 

where A, B, C, A are complex numbers, let us describe the asymptotic ratio ^f(z) = 
linin^oo p " + /f f ^ and the asymptotic ratio measure v for i?/- and i?//-type polyno- 
mial families. Using our previous results we can restrict ourselves to considering a 
recurrence relation of the form 

P„+l(z) = Ql{z)p n {z) + Q2{z)p n ~l{z), (4.1) 

where p-i(z) = 0, po(z) =1, degQi(z) = 1 and degQ2(z) < 2. For our study there 
is no essential difference between the two types. By Theorem[6]and Corollary [2] the 
asymptotic ratio *&(z) satisfies the symbol equation 

* 2 (z) = Qi{z)9(z) + Q 2 (z) 

everywhere in the complement of the union of the induced maxmod-discriminant 
E!$ and a finite set of (isolated) points. The support of v coincides with 2$. In 
order to describe S$ we need the following simple lemma. 

Lemma 5. A quadratic polynomial t 2 +ait+a2 has two roots with the same absolute 
value if and only if there exists a real number eg [1, oo) such that ea\ — 4a2 = 0. 



Proof. The roots of the polynomial are given by = ai± V a i 4a2 ^ g0 
that if |fi| — \t2 \ then the complex numbers a\ and ^Ja\ — 4a2 considered as vectors 
in C must be orthogonal. Thus | Arg(ai) — KvgyJ a\ — 4a2| = which is equivalent 

to saying that 4 ° 2 is purely imaginary and therefore ai ~f° 2 € (— oo, 0]. The 
converse statement is obvious. 

Corollary 4. The induced maxmod-discriminant 2$ C C is the set of all solutions 
z to the equation 

eQl(z) + AQ 2 (z) =0, (4.2) 

where e € [1, oo). 
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Figure 2 . Zeros of P41 (z) satisfying the 3-term recurrence relation 
p„+i(z) = zp n (z) + i^f^-)Pn-i(z) with p-i(z) = and p (z) = 1 

Below we concentrate on the situation when Qi(z) and Q 2 (z) are real polynomi- 
als. Let C denote the unique (real) root of Qi(z) and denote by D, E the branching 
points of the symbol equation, i.e., the roots of Q\{z) + 4Q 2 (-z) = 0. If Qi(z) and 
Q2{z) are real there are three basic cases: 

(1) D and E are real, D < E and C G [D, E]; 

(2) D and E are complex conjugate; 

(3) D and E are real, D < E and C <£ [D, E]. 

Case 1. The discriminant 5$ is the interval [D,E]. The density of v — in 
this case is positive and equals 



^Ql(x)+AQ 2 (x)dx, (4.3) 



where the value of yj Q 2 (x) + 4:Q2(x) is taken with negative imaginary part. (The 
value of Q\{x) + AQ 2 (x) on the interval [D, E] is negative.) 



Remark 7. Note that J D dv is not necessarily equal to 1. In fact, J D dv — 



(E-D) 
16 



For comparison let us take the sequence of inverse ratios ^—fe f° r the above 
sequence of monic polynomials {p n (z)}. Let v n denote the residue distribution 
of and set v = liirin^oo v n . One can easily see that the distribution v is 

actually a probability measure (i.e., it satisfies J~ dv = 1) having the same support 



as v and with density given by 



^Ql(x)+AQ 2 (x} dx 



2tt Q 2 {x) 

For the concrete example depicted in Figure [2] the asymptotic ratio measure v 
equals 2 ^~ X 2^^ ^dx and is supported on the interval [—4, 1]. 

Case 2. The discriminant S$ is the arc D — C — E of the unique circle passing 
through all three points D, C, E. In order to simplify the notation while calculating 
v let us assume that degQi(z) = 1. Then by an affine change of the variable z we 
can normalize Q\(z) such that Q\(z) = z. Set Q 2 (z) = az 2 + bz + c, where a, 6, c are 
real numbers. The condition that Q\{z) +4Q2(^) = has complex conjugate roots 
gives b 2 — c(4a + 1) < 0. Then the center of the circle passing through D, C, E is 
given by (— f , 0) and its radius is | f | . If 7 is the angle from the real positive half-axis 
to the ray through the center of the circle and a point (x, y) on it then we have the 
following parametrizations: x — ||| cos 7 — | , y = || | sin 7 and n = cos 7 + isin7 
for the unit normal to the circle at the point (x,y). 

In order to calculate the asymptotic ratio measure v = d gi 2 ^ recall the following 
statement from [10[ Theorem 3.19]. 

Proposition 3. Let Y C X be open subsets o/M fc such that Y has a C 1 -boundary 
dY in X and let u £ ^(X). If Xy denotes the characteristic function of Y , dS 
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Figure 3. Zeros of P41 (z) satisfying the 3-term recurrence relation 
p n +i(z) = zp n (z) + (i^)p„_i(z) with p-i(z) = and p (z) = 1 



the Euclidean surface measure on dY and n the interior unit normal to dY then 

dj( u X Y ) = i d j u )x Y +urijdS, 
where dj denotes the partial derivative with respect to the j-th coordinate and rij 
denotes the j-th component of n. 

From Proposition [3] we deduce the following result. 

Corollary 5. In the notation of Proposition^ one has 



dz 



du 



+ -u(ni + in 2 )ds. 



(4.4) 



Proof. Take an open set Y C M 2 with a C 1 -boundary dY oriented in such 
a way that Y lies to the left of dY. If s denotes the natural (length) parameter 
and dY is parametrized by (x(s),y(s)) then (ff,ff) is the tangent unit vector 
and (— 4^, 4^) is the interior normal unit vector to dY at the point (x(s),y(s)). 



It follows that m + in 2 = + i% = i{% + iff), so that formula 



is a 



consequence of Proposition [3] and the fact that 



dz 2 I dx ^ 1 dy ) ' 



Corollary [5] implies in particular that if *&(z) is a piecewise analytic function with 
smooth curves separating the domains where ^S(z) coincides with a given analytic 



function then v 



92 



is concentrated on these separation curves. Additionally, 



at each smooth point where two branches \&i and ^ meet v satisfies the condition 



(with an appropriate choice of coorientation). Applying now this 



(^1-^2) 

result to the case under consideration we see that formula (|4.4p reads 



= d^jz) = y/Ql(z)+4Q 2 (z)n 
V dz 2ir 
Explicit computations (using Mathematica ) give 



V2, 



v = —j— (cos 7 + i sin 7) 5 yj c(2b 2 — (1 + 4a)c(l — cos -yfjd'y. 
For the concrete example considered in Figure [3] one has 
v = \Jl cos 7 (cos 7 + % sin 7) 5 c/7, 



(4.5) 



where 7 is the angle from the real positive half-axis to the ray emanating from 
( — f ; 0) = (1,0) (i.e., the center of the circle of radius 1 1 = 1 containing the points 
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D,C, E in this case) and passing through a variable point (x, y) lying on the left 
half-circle. 

Case 3. The discriminant S$ is the union of the interval [D, E] and the circle given 
(as in Case 2) by the equation x[x — xq) + y 2 = where x (respectively, y) is the 



real (respectively, imaginary) part of z and xq 



2c 
b 



The density of v = 



d*S!(z) 



is 



given by formula (|4.3j) on the interval and by formula (|4.5p on the circle. A concrete 
example is depicted in Figure 2] below. 
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Figure 4. Zeros of pn (z) satisfying the 3-term recurrence relation 
p„+i(z) = zp n {z) + (^±^)p n -i{z) with p-i(z) = and p (z) = 1 

If cither of the polynomials Qi(z) and Qi{z) has complex coefficients it seems 
difficult to give a more precise description of the support of v than the one obtained 
in Corollary |U Figure O below illustrates some possible forms of this support. 
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-2 -1.5 -1 -0.5 



D.25 -0.2 -0.15 -0.1 -0.05 



Figure 5. Zeros of pa(z) for the recurrence relations p n+ i(z) = 
zp n (z) + (Iz 2 +2z-l- 3I)p n -i(z) (left) and p n +i{z) = zp n (z) + 
((2-i> 2 + z + 4-7>„_i(z) (right) withp_i(z) = andpo(-z) = 1 

Note that in all examples in this section we took the standard initial polynomials 
P-i(z) = and po(z) — 1. Computer experiments show no isolated zeros of p n (z) in 
these cases. The following lemma shows that in the specific situation when p_i(z) = 
and po(z) = 1 the set of slowly growing initial conditions (cf. Definition [7]) is 
actually empty. 

Lemma 6. For any recurrence relation of the form (|4.1|) with initial values Pi(z) — 
and Po(z) — 1 the set E$ of slowly growing initial conditions is empty, comp. [4]. 

Proof. A recurrence relation of the form (|4.ip may be rewritten as 

(Pn+x{z)\ _ rpt \ ( Pn(z) 



Pn{z) 



T{z) 



Pn-l{z)J ' 



where T 



{z) = (Qf) Qf)\ 
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It is easily checked that for any z£C the vector (1,0)' € C 2 is not an eigenvector 
of T(z), which immediately implies the desired result. 

We should point out that as soon as one changes the initial pair (p_i(z),po( z )) 
the polynomials p n (z) immediately acquire additional isolated zeros away from 5$. 

5. Related topics and open problems 

5.1. Does Theorem Etu) hold in the algebraic and smooth categories, i.e. when 
the varying coefficients are algebraic functions of bounded degree or smooth? 

5.2. For finite recurrence relations with fixed polynomial coefficients and arbitrar- 
ily prescribed initial values the following fact observed in computer experiments 
needs a mathematical ground and a conceptual explanation. 

Conjecture. The zeros of any two consecutive polynomials p n+ \{z) and p n (z) in- 
terlace along the curve S$ for all sufficiently high degrees n. 

An example illustrating the interlacing phenomenon conjectured above is shown 
in Figured The large fat points in this picture are (as usual) the branching points. 
Normal size fat points are the zeros of pn(z) while thin points represent the zeros 
of p4o(z). One notices a typical interlacing pattern between these two sets of zeros. 



"0 




, , , , , , t 

-2 -1.5 -1 -0.5 0.5 1 1.5 

Figure 6. Zeros of Paq(z) and pn{z) satisfying the 4-term recur- 
rence relation p n +x (z) = (z+1 — I)p n (z) + (z + l)(z — I)p n -i(z) + 
(z 3 + 10)p n - 2 (z) with p (z) = z 6 - z 4 + I , p 1 (z) = z - I + 2 and 
p 2 {z) = (2 + 7> 2 -8 

Similar interlacing properties for complex zeros along specified curves were stud- 
ied in e.g. [2]. Some caution is actually required when defining the interlacing 
property since the zeros of p n {z) do not lie exactly on S$. In the case when S$ is 
a smooth curve we may proceed as follows. Identify some sufficiently small neigh- 
borhood ./V(S$) C of S$ with the normal bundle to S$ by equipping 7V(S$) 
with a projection onto S$ along the fibers which are small curvilinear segments 
transversal to Then one can say that two sets of points in A^(S$) interlace 
if their projections on S$ interlace in the usual sense. If S$ has singularities one 
should first remove some sufficiently small neighborhoods of these singularities and 
proceed in the above way on the remaining part of The above conjecture then 
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says that for any sufficiently small neighborhood 7V(S$) of S$ and any of its iden- 
tifications with the normal bundle to 2$ there exists no such that the interlacing 
property for the roots oi p n (z) and p n+ i(z) holds for all n > uq. 

5.3. ft is natural to ask what topological information about the induced maxmod- 
discriminant may be obtained from the coefficients of a given recurrence relation. 

Problem 1. Describe the topological properties of S$ depending on the coeffi- 
cients of the recurrence relation (say assumed fixed and algebraic) and establish, in 
particular, necessary and sufficient conditions for the set S$ to be compact. Are 
there any characteristic numbers associated with the singularities on 5$? 

5.4. Proposition [2] determines the support of the asymptotic ratio distribution v 
associated to a recurrence relation of the form (|1.7[) with fixed coefficients (i.e., 
independent of n). However, the actual support of v in the general case is still 
unknown. 

Problem 2. Determine the support of v in the case of varying coefficients. What is 
the analog of Lemma [3] (iv) for varying coefficients? Note that for fixed coefficients 
and generic initial data the divergence of the ratio actually takes place on an open 
everywhere dense subset of S$. Under the same requirements the ratio converges 
on an open dense subset of but to a "wrong" limit (in the sense that this limit 
is different from the leading spectral number). 

5.5. Paper [1 contains some intriguing combinatorial and geometric facts related 
to a discriminant surface which is similar to the maxmod-hypersurface intro- 
duced in Definition If 21 More precisely, the discriminant considered in op. cit. con- 
sists of polynomials with at least two zeros having the same real part. This dis- 
criminant is also related to quadratic differentials. 

Problem 3. Extend the results of pQ to the maxmod-hypersurface 

6. Appendix. On topology and geometry of maxmod-discriminants 

Definition 12. Denote by C Polk = {t k +a\t + . . .+ak} the set of all degree 
k monic polynomials with complex coefficients having at least two roots with the 
same absolute value. 

It is obvious from Definitions [U and [T^] that contains 3/. . 

Proposition 4. is a real semialgebraic hypersurface in Polk of degree at most 
(4fc — l)(4fc — 2) considered as a 2k -dimensional real linear space with the real and 
imaginary parts of all ai 's chosen as coordinates. Furthermore, the hypersurface 
Sfc is quasihomogeneous with quasihomogeneous weight equal to i for both the real 
and the imaginary parts of <n, i — 1, 2, . . . , k. 

Proof. In order to calculate the degree of Efc let us describe an algorithm 
giving the equation for the analytic continuation of This algorithm may be pre- 
sented as a superposition of two resultants. Let u = x+iy £ C and consider first the 
resultant of P(t) = t k + a x t k - 1 + ... + a k and P(ut) = {ut) k + ai(ut) k ~ 1 + ... + a k , 
which we denote by i?(u, a±, . . . , <!&). (Recall that all aj's are complex.) Clearly, 
R(u, ai, . . . , dfc) is a polynomial in the variable u. For any fixed value of u the resul- 
tant R{u, a±, . . . , a*;) vanishes if and only if P(t) and P(ut) (considered as polyno- 
mials in t) have a common zero. One can easily see that R(u, a±, . . . , a^) is divisible 
by (u — l) k and that the quotient R(u, a\, . . . , a/-) = R(u, a%, ... , &&)/ (u — l) k is 
coprime with u — 1. We want to find an equation for the set of all monic complex 
polynomials P(t) such that P(ut) and Pit) have a common zero for some « ^ 1 
with \u\ = 1. The classical rational parametrization of the unit circle is given by 
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u = ^ 1+ 'g2 , where 6 G M. The result of the substitution u — x + iy — m 

i?(it, ai, . . . , afe) is a complex- valued rational function of the real variable whose 
denominator equals to (1 + 9 2 ) 2k . By taking now the resultant of the real and 
imaginary parts of the numerator of the latter rational function one gets the re- 
quired algebraic equation for the analytic continuation of S/.. This recipe allows 
us also to calculate (an upper bound for) the degree of S&. It is easy to see that 
R(u, a\, . . . , a/-) is a polynomial of degree 2k in u and of degree at most 2k — 1 in 

the variables a\, . . . , a k - After making the substitution u = x + iy = ^ 1+ 'gj and 
taking the real and imaginary parts of the above rational function one gets two 
polynomials of degrees at most 2k — 1 in the variables 5iai, . . . , SRa^, Sai, . . . , Sa^ 
and of degrees 4fc and 4fc — 1, respectively, in 6. These polynomials have propor- 
tional leading terms and their resultant is therefore a polynomial in the variables 
Rai,...,Ua k ;^ai,...,^a k of degree at most (2(4fc)-2)(2fc-l) = (4fc- l)(4fc-2). 

The quasihomogeneity of with quasihomogeneous weights as specified in the 
statement of the proposition follows from the fact that 5^ is preserved if one mul- 
tiplies all roots of a polynomial by a nonnegative real number. 

Remark 8. The proof of Proposition U contains an algorithm giving the explicit 
equation for the analytic continuation of Sj, . The authors wrote a Mathematica™ 
code doing this for small values of k. However, the resulting expression contains 
several hundred terms even for k = 3 and does not seem to be of much use. For 
k = 3 the estimate (4fc — l)(4fe — 2) for the degree is sharp. It is very plausible that 
this estimate is actually sharp for arbitrary k although in the general case it seems 
very difficult to check the necessary nondegeneracy conditions while performing the 
superposition of the two resultants described in the above algorithm. 

Proposition 5. Ejj. is a singular real hypersurface with a (singular) boundary coin- 
ciding with the usual discriminant T>k C Polk , that is, the set of all polynomials in 
Polk with multiple roots. Moreover, 2^ has no local singularities (i.e., singularities 
on a given branch) outside T>k- The list of singularities of 5^ is finite for any given 
k G N, i.e. its singularities have no moduli. 

Proof. Indeed, to show that T>k C Polk is the boundary of S/. consider the 
standard Vieta map Vi : C fe — > Polk sending a fe-tuple of (labeled) roots to the coef- 
ficients of the monic polynomial with these roots, i.e, to the elementary symmetric 
functions of these roots with alternating signs. It is known that the Vieta map 
induces a local diffeomorphism on the complement C k \7k, where Tk is the stan- 
dard Coxeter hyperplane arrangement consisting of ( 2 ) hyperplanes L^j given by 
Xi = Xj, 1 < i < j < k. An easy observation is that Tk coincides with the preimage 
V r i~ 1 (I?fe) of the discriminant T>k- Consider now the arrangement of quadratic cones 
£fc = Ui <:) C;.j in C , where is given by the equation \xi\ = \xj\. Obviously, Cfe 
coincides with the preimage Vi~~ l {z>k)- From the defining equation it is clear that 
each £ij is smooth outside the origin. Therefore, 5^ is locally smooth (that is, it 
consists of smooth local branches) outside T>k ■ Note also that for 1 < i < j < k the 
quadratic cone £jj contains the complex hyperplane Lij and that the restriction 
Vi\ € is a fold near a generic point of Lij. Therefore, T>k is the boundary of 

Efe . The absence of moduli for the singularities of can be derived from that for 
the singularities of T>k- The type of a (multi)singularity of near a polynomial 
P(z) € Efc is encoded in the following combinatorial information about the roots 
of P(z), which implies that the number of possible types of multisingularitics for 
a given k is finite. Determine first the multiplicity (possibly vanishing) of the root 
of P(z) at the origin and the set of all distinct positive absolute values for all the 
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roots of P(z). For each such positive absolute value determine how many distinct 
roots have exactly this absolute value and what are their multiplicities. If all the 
roots are simple then (as above) locally consists of (in general, nontransversal) 
smooth branches. Otherwise P(z) lies in T>k, which is the boundary of Efc, and 
local branches of near P{z) might have boundary singularities. 

Let be an open subset of M. 2 and assume that the coefficients of the recurrence 
relation (|1.7[) are sufficiently generic. The arguments in the proof of Proposition [5] 
applied to the induced maxmod-discriminant 5$ of the associated symbol equa- 
tion (|1.8|> imply that the only possible singularities on the curve S$ are the end- 
points and "Y-type" singularities, i.e., triple (local) rays emanating from a given 
point, see Figures HHS1 Note however that the transversal intersection of smooth 
branches (that is, the interval and the circle) on Figure[5]is actually unstable since it 
disappears under a small perturbation of the coefficients of the recurrence relation. 

Notation 4. Given a topological space X let X be its one-point compactification. 

Proposition 6. The following properties hold for any k £ N : 

(i) The hypersurfaces and are contractible in Polk, so that their usual 
(co)homology groups are trivial. 

(ii) The one-point compactification is Alexander dual in the sphere S 2k to a 
circle S 1 . Therefore, Sfc is homotopically equivalent to a sphere S 2k ~ 2 , so 
that i/jfSfcjZJ = Z if i S {0, 2k — 2} and iJj(Sfe,ZJ = otherwise. 

(iii) The one-point compactification is Alexander dual in the sphere S 2k to 
a (k- I) -dimensional torus T k ~ x . Hence Hi(S k) Z) ~ ff 2fe-i_1 (T fc_1 ,Z) 

for any i. Here Hi(X) (respectively, H l (X) ) stands for reduced mod point 
homology (respectively, cohomology) of the topological space X , see e.g. [26j . 

Proof. The contractibility of and Sj, follows directly from their quasiho- 
mogeneity. In order to show (ii) and (iii) we use the standard Alexander duality in 
Polk- Let Tfe = Polk \ Sfe and Tk = Polk \ Sfc denote the complements in Polk of 
Sfc and respectively. Then 

H 2k - 1 (T fe , Z) = Hi (S k , Z) and H 2k - 1 (T k , Z) = H. t (§ fc , Z) . 

Obviously, the space Ti ~ Ti ~ C is contractible. The next lemma describes the 
topology of Tfc and for k > 1. 

Lemma 7. For any k € N, k > 2, one has 

(i) Tk is an open 2k -dimensional manifold which is homotopically equivalent 
to the (k — 1)- dimensional torus T k ~ 1 . 

(ii) Tk is an open 2k -dimensional manifold Tk which is homotopically equiva- 
lent to a circle S . 

Proof, (i) The space Tk consists of all fc-tuples of complex numbers with 
distinct absolute values. Let Xk — {n < r 2 < ■ ■ ■ < rk \ r\ > 0} denote the set of 
all possible fc-tuples of distinct absolute values. Then is "fibered" over Xk with 
a "fiber" which is isomorphic to T k if r\ ^ and isomorphic to 7" fe_1 if r\ = 0. 
In order to get an actual fibration consider the set Xk — {0 < r% < . . . < rk} 
of the absolute values of the roots starting from the second smallest. Now T& is 
actually fibered over Xk with a fiber isomorphic to T fc_1 x D r2 , where D r2 stands 
for the open disk of radius r 2 > centered at the origin. The observation that Xk 
is contractible now implies that is homotopically equivalent to T fe_1 . 
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(ii) The space consists of all fc-tuples of complex numbers such there exists a 
unique number with largest absolute value in the considered fc-tuple. Let < r max 
denote this largest absolute value. Then Tfe is fibered over R + ~ {r max } with a 
fiber given by the product S 1 x Polk-x{r max ), where Polk-i{r max ) stands for the 
set of all polynomials of degree k — 1 whose roots lie in the open disk of radius r max 
centered at the origin. Since both R + and Polk-i{r max ) are contractible it follows 
that the space is homotopically equivalent to S 1 . 

Lemma [7] and the fact that S 1 is always unknotted in Polk for k > 2 complete 
the proof of Proposition [6] 
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